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Abstract 

The Jacobi equation in pseudo-Riemannian geometry determines 
the linearized geodesic flow. The linearization ignores the relative 
velocity of the geodesies. The generalized Jacobi equation takes the 
relative velocity into account; that is, when the geodesies are neighbor- 
ing but their relative velocity is arbitrary the corresponding geodesic 
deviation equation is the generalized Jacobi equation. The Hamilto- 
nian structure of this nonlinear equation is analyzed in this paper. 
The tidal accelerations for test particles in the field of a plane gravita- 
tional wave and the exterior field of a rotating mass are investigated. 
In the latter case, the existence of an attractor of uniform relative ra- 
dial motion with speed 2 -1 / 2 c ~ 0.7c is pointed out. The astrophysical 
implication of this result for the terminal speed of a relativistic jet is 
briefly explored. 

PACS numbers: 04.20.Cv, 98.58.Fd; Keywords: general relativity, Jacobi 
equation, jets 
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1 Introduction 



The analysis of observations in Einstein's general relativity theory is a rather 
complicated issue in contrast with the elegant simplicity of the geometric 
structure of the theory 0. Take, for instance, the excess perihelion preces- 
sion of Mercury that was explained by Einstein in 1915 and provided the first 
major success of general relativity. This standard theoretical result follows 
from the solution of the geodesic equation in the post-Newtonian approxi- 
mation. However, the relevant observational results are usually obtained by 
monitoring the motion of Mercury from the Earth via electromagnetic ra- 
diation that is reflected by Mercury and reaches astronomical telescopes on 
the Earth. In effect, one has to deal with the relative motion of one geodesic 
("Mercury") with respect to the other ("Earth") 0. Ideally, the observer 
following the reference geodesic sets up in its neighborhood a laboratory, 
where the measuring devices are assumed to work as in inertial spacetime. 
This local inertial frame of the observer is represented by the Fermi coordi- 
nate system. To construct such a coordinate patch, an orthonormal set of 
ideal gyroscopes would be required to define the local spatial frame of the 
observer, while an ideal standard clock would measure the observer's proper 
time. The reduced equation of geodesic motion of a nearby particle in such a 
system is the generalized Jacobi equation that will be studied in this paper. 
The results may be applied to high energy astrophysical phenomena as well 
as precision measurements in space involving the tidal field of the Earth or 
an incident gravitational wave. 

In a congruence of timelike geodesies, neighboring curves have rates of 
separation that are usually negligible compared to the speed of light in vac- 
uum c. In this case the geodesic deviation equation is the standard Jacobi 
equation |Q. There are physical circumstances, however, where the worldlines 
of adjacent geodesies diverge rapidly. Consider, for instance, a congruence 
of test particles with a common initial velocity falling toward a black hole 
such that depending on their impact parameters, some particles would be 
captured by the black hole, while others would escape its gravitational pull. 
Though the geodesies are neighboring as they approach the black hole, the 
rate of separation between the particles that fall in and those that are merely 
deflected would no longer be small compared to c. This example makes it 
clear that the relevant equation of relative motion, i.e. the generalized Ja- 
cobi equation, has a limited domain of applicability, since rapidly diverging 
geodesies would no longer be adjacent after a period of time that is char- 
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acteristic of the spacetime curvature — that is, g/c, where g is the radius of 
curvature. 

The investigation of the relative motion of neighboring observers follow- 
ing arbitrary timelike curves leads to the general deviation equation. The 
relativistic theory of tides is based on the general geodesic deviation equa- 
tion. This is rather complicated when the rate of separation of geodesies 
could be close to the speed of light; therefore, we confine our detailed discus- 
sion here to a limited form of this equation that is valid to first order in the 
relative separation (i.e. the generalized Jacobi equation). On the other hand, 
in the ( "nonrelativistic" ) case where a smooth geodesic congruence is given, 
a general treatment of the geodesic deviation equation is available; in fact, 
a new derivation can be found in the recent work of Kerner, van Holten and 
Colistete f§] who also give a novel application: they use the geodesic devia- 
tion equation together with Poincare's perturbation method to approximate 
bounded (eccentric) geodesies in the Schwarzschild metric that are close to 
circular geodesies. 

The generalized Jacobi equation is discussed in detail in Section 2; in 
particular, we are interested in the dynamical system represented by this 
nonlinear equation, where the nonlinearity is due to the existence of velocity- 
dependent terms in the system. Therefore, we present a detailed analysis of 
its Hamiltonian structure. 

To illustrate certain general dynamical features of the generalized Jacobi 
equation, we consider the relative motion of neighboring test particles in the 
field of a plane gravitational wave in Section 3 and in the exterior field of a 
gravitating source in Section 4. We note that many compact astrophysical 
systems emit powerful oppositely-directed jets that are highly collimated. 
No basic theory of astrophysical jets is available at present, though many — 
especially magnetohydrodynamic — models have been proposed. In particu- 
lar, the physical processes that determine the asymptotic speed of the bulk 
flow are not clear at this time. Of particular interest in this connection is 
the relevance of the relative radial motion described in Section 4 to the "ter- 
minal" bulk motion of a relativistic jet once such an outflow is beyond the 
active region around the source and the bulk motion is subject only to the 
gravitational attraction of the source. Taking only the relativistic tidal forces 
into account, we demonstrate in this case the presence of an "attractor" at a 
Lorentz factor T = y/2, corresponding to V — w 0.7c, for the terminal 

bulk motion of a jet relative to the ambient medium. 

Finally, Section 5 contains a discussion of our results. In the following 
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sections, we use units such that c = 1. 



2 Generalized Jacobi equation 

The generalized Jacobi equation in Fermi coordinates can be simply obtained 
by writing an alternative form of the geodesic equation, i.e. the reduced 
geodesic equation ( |2.5| ) below, in Fermi coordinates. However, to clarify the 
Hamiltonian and Lagrangian aspects of this equation, it is useful to proceed 
as in the following subsections. 



2.1 Hamiltonian and Lagrangian formalisms 

As is well-known, the geodesic equation in local coordinates x M on a pseudo- 
Riemannian manifold with metric 

ds 2 = g^(x)dx fi dx" 

is equivalent to the Euler-Lagrange equation for the Lagrangian 

£(/,/) = ^ V (^V. 

where 

dx^ 

U ■ dA' 
the canonical momentum is given by 

and A is an arbitrary parameter (see P). In fact, the Euler-Lagrange equa- 
tion for this system is equivalent to 

d 2 x fJl dx a dx 13 

d\ 2 a ? dX dX ' 

where 
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are the Christoffel symbols. The corresponding Hamiltonian is 

which is independent of the parameter A; thus, H is constant on geodesies. 
It follows that along a geodesic A is a linear function of s. 



2.2 Isoenergetic reduction 

Let us specialize to a four-dimensional Lorentzian manifold where we take 
an admissible coordinate chart and the sign convention (— , +, +, +) for the 
Lorentzian metric. We will also use the Einstein summation rules; Greek in- 
dices run from zero to three and Latin indices run from one to three. More- 
over, round brackets around indices denote symmetrization, while square 
brackets denote antisymmetrization. 

The timelike geodesies ((ds/d\) 2 < 0) with energy H — —1/2 have A = t, 
where r is the proper time along the geodesic such that (ds/dr) 2 = — 1. 
We will determine the isoenergetic reduction for the Hamiltonian system 
restricted to this manifold. 

The Hamiltonian equations of motion (where we have assumed for sim- 
plicity that p^ = w M ; that is, the mass of the free test particle is normalized 
to unity) are given by 

dx^ 1 &H „ 
dr dpf, 

d PlI OH _ 1 p 

~d^ ~ ~d^~~2 9 (2 ' 1} 

Reduction is accomplished by first separating the pairs of canonically con- 
jugate variables (x^,p v ) into time and space variables; that is, we take 
(x^,p u ) = (t,p ;x\ P j) so that points in spacetime with energy —1/2 sat- 
isfy the equation 

H(t,p 0] x\ Pj ) = ^g 00 p 2 + g 0l p oPt + ^g ij PiPj = ~. (2.2) 

Because the flow of time does not stop along the worldline of an observer, 
we have dt/dr ^ 0. Hence, by Hamilton's equations, the partial derivative 
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dH I dpo does not vanish. By the implicit function theorem, there is a function 
a such that 

H{t,a(t,x i ,p j ),x i ,p j ) = ~. (2.3) 

That is, the energy surface is (locally) the graph of a. Of course, we can also 
solve explicitly for a. Indeed, because g 00 ^ in an admissible chart we can 
complete the square with respect to p in equation ( |2.2|) to obtain 

,0i 



where 



Hence, 



\p° + ~^m> - ^oo 



qV ■- q v _ 



a(t,x, Pj ) = —^ Pi ±{ ^_ g0Q) ) 



We note that (g 4J ) is the inverse of the spatial metric (<7y); that is, <7* fe y";y = 5*. 
On the energy surface, the equation 

dpo dU f , u i \ i \ 
— = --Qjr(t,a{t,x ,pj),x ,p d ) 

decouples from the system ( |2.1| ). By computing the partial derivatives with 
respect to x % and Pi on both sides of equation Q2.3Q , we obtain the identities 

dH da_ dH 

dpo dx l dx l 
&Hdc^ + dH _ Q 
dp dpi dpi 



Using these relations, the Hamiltonian system (271) for the spatial variables 
can be expressed as 

dx l dH dH da 



dr dpi dp dpi 

dpi dH dH da 

dr dx l dpo dx l 
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where dTi/dpo = dt/dr 7^ since t is the timelike variable. Hence, the 
remaining "spatial" equations in the Hamiltonian system fl2.1| ) can be recast 
into the time-dependent Hamiltonian system ("the isoenergetic reduction") 

dx 1 da a . 

— = - Wi (t,x', Pk ), 

f = f^.ft). (2.4) 

with Hamiltonian —a. 

The Hamiltonian isoenergetic reduction procedure is equivalent to a La- 
grangian procedure. In fact, let us define the Lagrangian 

L(t, x\ v j ) = -[- (g 00 (t, x) + 2g 0i (t, x)v* + g^t, x)«V)] 1/2 , 

where v % := dx l /dt, the minus sign under the square root ensures that L is 
real and the other minus sign is to conform with the sign convention used in 
the Hamiltonian system ( |2.4| ). We now have 

dL 

Pi-= q~ = (m + QijV 3 ) [ - (#oo + ZgoiV 1 + 9ijV l v nl ' ~ 
and the corresponding time-dependent Hamiltonian 

H = H(t, x\pj) -.^PiV 1 - L(t,x\v j ) 



is given by 



TT r doc ^ dt 

H = -{g o + 9Qi-rr)-r 
dt dr 



= -(goo? + goiP 1 ) = -goaP a = -po, 

where we have used the fact that the proper time r along the geodesic world- 
line is such that L = —dr/dt. Hence H = —a in system ( p.4| ). 

This reduced system is equivalent to the second-order differential equation 

d 2 x l . dx a dx 13 .dx 1 t dx a dx 13 . . 

~w ~ { ^~df~df } ~dt + a(3 ^r^r = °' (2 - 5) 

which is the reduced geodesic equation. To prove this fact, note that the 
Hamiltonian or Lagrangian equations of motion are equivalent to the Euler- 
Lagrange equation 

j t [(goi + 9ijV j )[-{goo + 2g 0i v i + s^V)] -1 / 2 ] 
= -{g m + 2g ojti v j + g k i,iV k v e )[-{g 00 + 2g 0i v i + g ijV l v j )]~ 1/2 . (2.6) 
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Let us use the proper time r to recast equation ( |2.6| ) in the form 

^(9i a u a ) = ^g a p,iU a u p , 

where = dx^/dr. After computation of the derivative by the product rule, 
the equation can be rearranged to the equivalent form 

du a 1 

9ia-r~ + 7i(9ia,l3 + 9if3,a ~ 9 'af3 ,i)U° l vP = 0. 

We multiply both sides of the last equation by g^ % and use the identity 
9^9ia + g^goa = K to show that 

du a 1 

(<$a - g^goa)-^- + ^(gia,P + 9ip,a - 9afi,i)U a vP = 0. 

Also, using the identity 



— ^g^igiafi + 9i/3,a ~ 9a/3,i) + l^^^Oa^ + #0/3,0 ~ 9ap,o), 



we have the equation 

— + v%u a vP - gP&Jjp - ^°(goa, P + go?,* - g aPfi )u a u p = o. (2.7) 

We claim that 

M^ + r> a ^) = o. (2.8) 

To prove this identity, we replace F^a by its defined value, multiply through 
by Ufj, and then use the antisymmetry of the last two terms in the resulting 
expression to conclude that 



.du* 1 „„ „, As du^ 1 1 



— + r>«V) = Ufl — + -u' w «V. 

A calculation shows that the right-hand side of this last equation is equal to 
and therefore it is equal to zero, since (ds/dr) 2 = g^u fJ, u u = —1. 
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After multiplying equation (|2.7| ) by u^, we use equation ([2.8[ ) to see that 
c/m° 1 

fi'oa-^r + 2^ 0q ' /3 + ~ 9 a (3,o)u a u p = 0, 

since = u° = dt/dr ^ 0. By inserting this equality into equation ( |2.7|) , 

it follows that 

dv^ 



In particular, if /x = 0, then 



^ + I>V = 0; (2.9) 



and if \x = i, then 

d 2 x % 



dr 2 ' ~ a/3 
Using the chain rule, we have that 



+ r i ap u a u> 3 = 0. (2.10) 



ciV / x 2 (i 2 x* d 2 t dx l 



dr 2 dr dt 2 dr 2 dt 



By inserting this identity into equation ( 2.10 ) and by replacing dx p jdr with 
(dx p /dt) (dt/dr), we obtain equation ( |2.5| ) with its left-hand side multiplied 
by the nonzero factor (dt/dr) 2 . This proves the desired result. 



2.3 Generalized Jacobi equation 

The generalized Jacobi equation is the linearization of the reduced geodesic 
equation ( |2.5| ) in Fermi coordinates with respect to the space variables, but 
not their velocities. We will obtain the explicit expression for this equation in 
Fermi coordinates valid in a cylindrical spacetime region along the worldline 
of a reference observer, which is at rest at the spatial origin of the Fermi 
coordinates (see Appendix A). This is perhaps the most useful form of the 
generalized Jacobi equation. 

For an arbitrary spacetime, when a reference observer (that follows a 
geodesic) in a coordinate chart (t, x % ) is specified, we can always find a Fermi 
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coordinate system (T, X*) in a neighborhood of the worldline of the fiducial 
observer so that the metric tensor is given by 

Soo = -l- F Ro iaj (T)X i X S + -'-> 
9oi = ~ F R ojik (T)X j X k + --- , 

9a = - \ F R im {T)X k X t + • • • , (2.11) 



where 



F Raf3~{& — ^^A^A^A^A^), 



are the components of the Riemann tensor projected onto the orthonor- 
mal tetrad frame A^(T) that is parallel propagated along the worldline 
of the reference observer at rest at the spatial origin of the Fermi chart (cf. 
Appendix A). By expanding the coefficients of the reduced geodesic equa- 
tion ( p.5|) for the metric expressed in the Fermi coordinates and retaining 
only the first-order terms in the spatial directions, we obtain the generalized 
Jacobi equation 



d 2 X l 



+ F RoiOjX j + 2 F R ik joV k X j 



dT 2 

+(2 F R Qm V i V k + F R ikje V k V e + 2 - F R Qkji V i V k V^ = 0, (2.12) 

where V 1 := dX % jdT . Alternatively, we could start from the Hamiltonian 
equations ( |2.4| ) and substitute in the Hamiltonian for g 00 , g° l and g l i expan- 
sions corresponding to equations ( |2.11| ) and recover equation ( |2.12| ) directly 



from Hamilton's equations ( |2.4| ); this is done in Appendix B. The generalized 
Jacobi equation in arbitrary spacetime coordinates is given in Appendix C. 



It is interesting to note a special scaling property of equation (|2.12|) . To 
this end, we let g denote the radius of curvature; it is defined here so that 
1/ q 2 is the supremum of the absolute magnitudes of the components of the 
Riemann tensor over the coordinate patch where the Fermi coordinates are 
defined. We will use the natural rescaling of the generalized Jacobi equation 
given by 

X 1 = qx\ T = gl (2.13) 
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In fact, under this rescaling the generalized Jacobi equation maintains its 
form except that each component of the Riemann tensor is multiplied by g 2 . 
The scaled equation is valid for |ir l | <C 1 and for arbitrary velocities whose 
magnitudes do not exceed the speed of light c — 1. The requirement that 
|x| <C 1 would naturally imply that in most cases the generalized Jacobi 
equation would hold for a limited interval of time t — to starting from some 
initial state at t = to such that \t — t \ < 1. 

The Hamiltonian dynamics of the generalized Jacobi equation will be 
explored in the following section for relative motion in the field of a plane 
gravitational wave. The dimensionless Fermi coordinates (t, x) introduced 
in equation ( 2.13 ) will only be employed in the latter part of Section 3 and 
should not be confused with the usual local coordinates that are used to 
express the spacetime metric of the plane wave below. 



3 Generalized Jacobi equation for a plane grav- 
itational wave 

In this section we will determine the generalized Jacobi equation for a plane 
gravitational wave given by the metric 

d s 2 = -dt 2 + dx 2 + U 2 (e 2h dy 2 + e~ 2h dz 2 ) } 

where U and h are functions of retarded time u = t — x. This metric is a 
solution of Einstein's field equations in vacuum (R^ = 0) provided that 

U, uu + (h u ) 2 U = 0. (3.14) 

Under this assumption, the gravitational field is of Petrov type N (see ||) 
and represents a plane wave propagating along the positive x-axis. 

When \h\ <C 1 and U — 1, the metric reduces to a linear plane gravita- 
tional wave in the T-T gauge that is linearly polarized ("©"). 

Any observer that remains at rest at spatial coordinates (x, y, z) follows 
a geodesic of this spacetime. This is a consequence of a more general result 
discussed in Appendix D. 

Let us consider the reference observer at rest at the spatial origin. We 
claim that the orthonormal tetrad frame of this observer 

Af 0) = (1,0,0,0), Af 1) = (0,1,0,0), 

Af 2 , = (0, 0, U- X e~ h , 0), Af 3) = (0, 0, 0, ET 1 ^) (3.15) 
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is parallel propagated along its worldline; that is, 



r 1 /vr^en^Crtl — u - 



dt P a (°) ( a ) 

It follows from Appendix D that the worldline of this observer lies on a 
geodesic (that is, the last equation holds for a = 0), so it suffices to show 
that the equation holds for a — i. In fact, because the nonzero components 
of T l 0a are 

F 2 1 #22,0 p 3 1 #33,0 

it is easy to see that A^ satisfies the equation. The desired result for and 
A^ follows from the fact that U~ x exp(— h) and U~ x exp(h) in equation ( |3.15| ) 

are respectively equal to (f?22)~ 1,/2 and (^s) -1 ^ 2 up to a sign factor. 

To obtain the generalized Jacobi equation for the reference observer, we 
must compute the Fermi components of curvature along its worldline. Since 
the observer is at rest at the origin (x = y = z = 0), these components will 
be functions of its proper time t = T only. In fact, for 

K(u) := h uu + 2h : J^j- 

we have that 

F R\212 — F Rl22Q — ^2020 = ~K(T), 
F ^1313 = ^1330 = ^3030 = K(T) 

are the only nonzero components except for the symmetries of the Riemann 
tensor. In the Fermi system, the spatial coordinates X are characterized 
invariantly with respect to the tetrad system as described in Appendix A. 

By inserting the curvature components into the generalized Jacobi equa- 
tion ( |2.12| ) with X' L = dX 1 / dT, we have the following system of second-order 
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differential equations: 



X = --K(T)(X - l)[X{Y 2 - Z 2 ) - YY{X - 3) 
3 

+ ZZ(X-3)}, 

Y = --K(T){2XY(X + Y 2 - Z 2 ) - Y[3 + 2Y 2 (X - 3) - 6X + 2X 2 ] 
3 

+ 2ZFZ(X-3)}, 

Z = -K(T){2XZ(X -Y 2 + Z 2 ) + 2YYZ(X -3) 
3 

-Z[3 + 2Z 2 (X-3)-6A > + 2X 2 ]}. (3.16) 

We note that this system is invariant under the transformation IhI, 
Y i— > Z, Z i— > y and X t— > —if. The generalized Jacobi system ( |3.16| ) with its 
velocity-dependent nonlinearities may conceivably have practical applications 
in the future for the case of space-based gravitational wave detectors such as 

LISA m. 



3.1 Dynamics 

System (|3.16| ) has a seven- dimensional extended state space such that the 
five-dimensional submanifolds 

M Y := {(X,Y,Z,X,Y,Z,T) :Y = Y = 0}, 
M z := {(X,Y,Z,X,Y,Z,T):Z = Z = 0} 

are invariant. The extended state space is simply the space of positions, 
velocities and time. On the three-dimensional invariant intersection of these 



manifolds, the flow is given by X = and T — 1. That is, equations (|3.16|) 
permit the uniform motion of a particle along the direction of wave propaga- 
tion relative to the fiducial particle. This is due to the transverse character 
of the gravitational wave. Since the corresponding two-parameter family 
of solutions is X(T) = X(0)T + X(0), the relative motion is always un- 
stable. That is, most members of this two-parameter family of geodesies 
diverge from the reference geodesic, which is the member of this family with 
X(0)=X(0) =0. 

In the small-/;, approximation for a monochromatic wave with U — 1, 
the curvature is given by K = h )UU ; hence, the assumption that h{u) = 
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— e cos((p+Qu) results in the curvature K(T) = eQ 2 cos(4>+QT) , where e< 1, 
Q is the wave frequency and is a constant phase. In this case, the radius of 
curvature is g = 1/ (£ly/e) and, using the scaled variables ( |2.13| ), system ( |3.16| ) 
maintains its form except that K(T) is replaced by k(t) := cos(0 + £///), 
where /i := sfi. It is convenient therefore to define a new temporal variable 
S :=t/fi and consider the equivalent first-order system with 

& = v£, y' = w, z' = fj-C, 

where a prime denotes differentiation with respect to the new "fast" time S. 
Each component of the vector field corresponding to the resulting first-order 
system has order /i and is 27r-periodic. This system is therefore in the correct 
form for averaging with respect to the temporal variable. 

By averaging to second-order in the small parameter /i, we transform 
system ( |3.16|) (by a near-identity nonautonomous change of variables) to a 
system that is autonomous to second-order in //. In fact, the second-order 
approximation of the averaged system is given by 

x' = fj£, 

£' = 0, 

y = /j,t] - fi 2 y sin 0, 

rj' = /j 2 ^ sin (f), 

z' = fi( + fi 2 z sin (f>, 

(' = -/i 2 Csin0. (3.17) 

This second-order approximation is not affected by the nonlinear terms in 
the generalized Jacobi equation; that is, the result would be the same for the 
corresponding Jacobi equation. The effect of the nonlinearities first appears 
in the third-order averaged system, i.e. at order e 3 / 2 . 

Let us note that taking into account the first-order approximation in 
[A only, the approximation procedure suggests that on average X, Y and 
Z are not affected by the presence of the wave; that is, to first order in 
fi, X(T) = X(0) + X(0)T. In case sin0 ^ 0, the first correction to this 
approximation — necessary only for Y(T) and Z(T) — is obtained from the 
second-order approximation of the averaged system (|3.17 ); in fact, we have 
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that 



Y(T) = Y(0)e 



-eQT sin < 



+ 



F(0) 



sinh(ef2T sin < 



ef2 sin < 

Y (0) + y (0)T - eF (O)fiT sin + 0(e 2 ) 



Z(T) = Z(0)e' 



e^Tsinc 



+ 



Z(0) 



sinh(ef2T sin . 



efi sin i 

Z(0) + Z(0)T + eZ(0)nT sin + 0(e 2 ). 



(3.18) 



It follows from these results that for sin0 = 0, the motion is on average 
unaffected by the presence of the wave to first order in e, as expected from 
the form of equation ( 3.17 ). Moreover, equations ( 3.18|) ignore nonlinearities 



in the system (|3.16| ) and are therefore valid for the Jacobi equation averaged 
to order e. It should be noted that for the Jacobi equation the motions in Y 
and Z are stable as follows from an analysis of the corresponding Mathieu 
equations. This is not in conflict with equations (|3.18j ) since to see this 
stability one has to average to higher orders than /i 2 = e. 



3.2 Hamiltonian Chaos 

The dynamical behavior of the generalized Jacobi system ( |3.16| ) is compli- 
cated by its high dimensionality and the velocity-dependent nonlinear terms. 
On the other hand, the system is Hamiltonian and if K is periodic, we ex- 
pect some dynamical effects associated with Hamiltonian chaos: resonance 
phenomena, stochastic regions and diffusion. While an investigation of these 
phenomena is beyond the scope of this paper, we illustrate the existence of 
Hamiltonian chaos in generalized Jacobi equations for the simpler system 
with 1 + 1/2 degrees of freedom given by 

d ^ = K(T)(l-2V*)X, (3.19) 

where V := dX/dT and K is a periodic function. This equation of rel- 
ative motion is obtained from the generalized Jacobi equation by limiting 
attention to a two-dimensional world, i.e. a surface (T, X), where K is the 
Gaussian curvature of the surface (cf. Section 4). According to the results 
of Appendix B, equation ( p,19|) can be obtained from a Hamiltonian 



H(T,X,p) = T (l +V 2 ) l/2 [l + \k{T)X 2 } 
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once only terms linear in X are retained in the equation of relative motion. 
It is interesting to remark that the equivalent first-order system 



X = V, V = K(T)(1-2V 2 )X, 



(3.20) 



is Hamiltonian with respect to the symplectic form 



l 
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Figure 1: The plot of V versus X at T = 2mi for n = 0, 1, 2, ... , 600 is 
depicted for several orbits of the Poincare map for the system ( 3.19Q with 
K(T) = 0.5cos(T). 



■\-2V 2 Y l/2 dX AdV, 



in the region of the state space where V 2 ^ 1/2. Also, we note that the 
change of coordinates 

for V 2 7^ 1/2, transforms system (|3.20|) to 

J_/ exp(2^)-l N . = _ K(T)x 
^Vexp(2v / 2^) + l^' ' 

a system that is Hamiltonian with respect to the usual symplectic structure. 
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A few typical orbits with \V\ < (1/2) 1 / 2 are illustrated in Figure |I| 
for a stroboscopic Poincare map obtained by numerical integration of sys- 
tem Q3.20D with K{T) = 0.5cos(T); that is, we plot V versus X at times 
2n7T for nonnegative integer values of n. We note the presence of stochas- 
tic layers, resonant island chains and invariant closed curves. The proof of 
the existence of chaotic invariant sets even in this two-dimensional case is 
not simple since a straightforward application of the Melnikov method runs 
into standard difficulties; in particular, one has to deal with the problem of 
exponentially small splittings of homoclinic loops. 

4 Generalized Jacobi equation for radial mo- 
tion in black hole spacetimes 

Imagine the relative motion of nearby free test particles along the axis of ro- 
tational symmetry of a Kerr source. For the reference geodesic, the equation 
of motion, using the standard Boyer-Lindquist radial coordinate r, is 0] 

where G is Newton's gravitational constant, M is the mass of the Kerr source 
and aM is its angular momentum. Here 7 is a positive constant of integration 
and 7 > 1 has the interpretation of the Lorentz factor for the particle at 
infinity (r — > 00). 

For a neighboring geodesic, the generalized Jacobi equation describing 
motion along the axis relative to the fiducial observer is 

^ + MT)[l-2(f)>^0, (4.22) 

where X = (R, 0,0), R is the relative radial distance along the axis, T = r 
along the reference geodesic and k(T) is obtained from the solution of 
equation ( |4.21| ) and 



_ GAfr(r»-3a» ; 



It is interesting to note that in this form k is independent of boosts along the 
axis. This is a consequence of the fact that the Kerr solution is of Petrov type 
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D and hence the Kerr symmetry axis provides two special tidal directions for 
ingoing and outgoing trajectories |8|, [|. For a = 0, the problem reduces to 
purely radial relative motion in the exterior Schwarzschild spacetime. 

More generally, starting from the generalized Jacobi equation and limiting 



our attention to a two-dimensional world as in Section |3.2| , let us consider 
the equation of relative motion in the form 

*L + K (i _ 2V 2 )X = 0, (4.24) 

where k(T), k = F i?oioi ; is the Gaussian curvature of the surface (X, T) and 
V = dX/dT. 



V 



x A 



X 



x_ 



Figure 2: Phase diagram of equation ( |4.24| ) for negative constant curvature. 
The bold horizontal lines correspond to X±. 

To investigate the main consequences of equation ( f4.24|) beyond the origi- 
nal Jacobi equation for arbitrary curvature, we assume that |V|, —1 < V < 1, 
is not negligibly small compared to unity. In particular, let us note that equa- 
tion (|4.24j) has constant velocity solutions with V 2 = 1/2 regardless of the 
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Figure 3: Phase diagram of equation Q4.24 ) for positive constant curvature. 
The bold horizontal lines correspond to X±. 



curvature of the surface, i.e. X± = X ± (1/2) l l 2 {T — T ), where T is a con- 
stant initial time and X = X(T ). It is therefore interesting to investigate 
the general behavior of solutions of equation ( |4.24|) near these constant ve- 
locity solutions. The asymptotic behavior of nearby solutions is determined 
by the curvature. Figures ^| and |3] depict the phase portraits of system ( [4.24]) 
for constant curvature cases k < and k > 0, respectively. We note that for 
constant k, equation ( |4.24| ) is completely integrable. It would be interesting 
to investigate the behavior of the solutions of equation ( f4.24|) near X± for 
arbitrary k(T). In the rest of this section we limit our attention to k(T) < 0, 
since it follows from equation ( |4.23|) that this is the relevant case for motion 
away from a Kerr source (r 2 > 3a 2 ). 

For the negative curvature case we consider the asymptotic stability of 
the invariant sets X± := {(X, V,T) : V = ±(1/2) 1//2 } in the extended state 
space for equation ( |4.24| ). For example, we call X + asymptotically stable if for 
every solution (X(T), V(T),T) in the extended state space of equation (|4.24|) 
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that starts near X + , its distance to X + along the velocity direction decreases 
to zero as T increases to infinity. It can be shown that asymptotic stability is 
essentially determined by the linearized first-order system along X + , i.e. the 
system involving the linear approximation W of X— X + and Q of V — (1/2) 1//2 
that is given by 

^ = Q, ^ = 2k(T)[V2X + (T-T )]Q. (4.25) 



By integration of system (|4.25|) , we find that the invariant set X + of the 



equation ( |4.24[ ) is positively asymptotically stable if for each real number Co 
we have 

/•oo 

/ k{T) (T + C ) dT = -oo. (4.26) 

This result has a useful corollary: if k(T) is asymptotic to a negative constant 
multiple of T~ p with p < 2 as T — ► oo, then X + is asymptotically stable. 
Let us note that if 7 < 1 in the differential equation ( 4.21 ), then there 



is a real solution only if the right-hand side of the differential equation is 
positive. In this case the real solutions with r > remain bounded between 
the two positive roots of the right-hand side. Under the additional assump- 
tion that the square of the smaller positive root exceeds 3a 2 , the curvature 
obtained from equation ( |4.23 ) is bounded below zero; therefore, the integral 



condition (|4.26| ) is satisfied and X + is positively asymptotically stable. For 
7 = 1, by using an asymptotic analysis of equation ( |L21j ) we find that r(T) 
is asymptotic to a constant multiple of T 2//3 as T — > 00 and k(T) is therefore 
asymptotic to a negative constant multiple of T~ 2 . Thus, for 7 = 1 the 
invariant set X + is again asymptotically stable; the rate of attraction, how- 
ever, is slower than any exponential rate. For 7 > 1, the set X + is no longer 
asymptotically stable. Although nearby solutions in this case move closer 
to X + as T — > 00, they are bounded away from this manifold. Let us note 
that the results of this analysis also apply to the behavior of solutions near 
X_, since in this case the linearized system is identical with equation (|4.25|) 
except for X \— > — X . 

The k < case is particularly interesting in connection with the problem 
of "superluminal" jets in astrophysics 0. In fact, the Galactic "superlumi- 
nal" jet sources GRS 1915+105 and GRO J1655— 40 have speeds that may 
be near (1/2) 1 / 2 w 0.7. 
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5 Discussion 



The main dynamical features of the generalized Jacobi equation have been 
presented in this paper. In contrast to the standard Jacobi equation, the 
generalized Jacobi equation is nonlinear and therefore exhibits Hamiltonian 
chaos under certain circumstances. Two specific applications have been con- 
sidered: relative test particle motion in the field of a plane gravitational wave 
and in the field of a rotating mass. The results may be relevant respectively 
to future space-based gravitational wave detectors and high-energy astro- 
physical phenomena associated with jets. 

The bulk speed of Galactic "superluminal" jets appears to be constant 
and estimated to be around 90% of the speed of light, though many uncer- 
tainties are associated with this estimate ||. Jet motion is detected via the 
temporal evolution of the radiation intensity contours relative to the back- 
ground, which would indicate the motion of blobs (electron clouds) in the jet 
relative to the ambient medium. This motion, under the influence of gravi- 
tation alone, is described by the generalized Jacobi equation as discussed in 
detail in the previous section. In this case, the generalized Jacobi equation 
exhibits an attractor with constant speed 2 -1 / 2 , approximately equal to 70% 
of the speed of light. 

Assuming that far from the source gravitational tidal forces are domi- 
nant, it follows from our results that the jet speed should approach 
from above and below over a certain characteristic timescale. This timescale 
can be calculated from the solution of equation Q4.25Q for Q. More explicitly, 
let us consider the jet along the axis of a Kerr black hole of mass M and 
angular momentum J < M 2 at a distance r ^> 2GM; then, the charac- 
teristic timescale is p , where p is the radius of curvature at r given by 
Po ~ (?"q/2GM) 1//2 . Simple estimates suggest that further refinements of the 
observational techniques B are required before it may be possible to test this 
prediction of the elementary theory of the asymptotic jet speed presented in 
this paper. 



A Fermi coordinate systems 

Consider a reference observer following a curve C in spacetime. Let r be the 
proper time along C and = dx^/dr be its tangent vector. The fiducial 
observer carries a triad A£\ representing three ideal gyroscope directions so 
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that A^ are parallel transported and A/\ with X^, = is an orthonormal 
tetrad frame. The Fermi system is a geodesic coordinate system valid in a 
cylindrical spacetime region around C such that for any point P with Fermi 
coordinates X a = (T, X % ) there exists a unique spacelike geodesic orthogonal 
to C that connects P with a point P on C. Let r be the proper time at 




- P 



Figure 4: Schematic diagram depicting the construction of the Fermi normal 
coordinate system. 

Pq, a be the proper length of the spacelike geodesic starting from P and 
rj^ = (dx^ I da) p be the vector tangent to the spacelike geodesic at Pq; then, 
f]^u^ = and the Fermi coordinates are defined by 

T = r, X^FarirXfa (A.l) 

(cf. Figure 3). 

To determine the metric tensor g^i^x) in Fermi coordinates, i.e. F g^ v (X), 
we expand F g^ v (P) in a Taylor series about F g^ v {Po) = f]^. In this expansion, 
the required derivatives of F g at P : (T, 0) are calculated on the basis of the 
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following considerations: The equations of motion of the tetrad frame carried 
by the observer along the fiducial geodesic C in Fermi coordinates reduces to 
f Tq (T, 0) = and F TQ i (T, 0) = 0. Moreover, the spacelike curve connecting 
Po to P satisfies the geodesic equation, which in Fermi coordinates (|A.l ) 
reduces to F T^ j (T,X)X i X^ = 0. Expanding F T^ in a Taylor series about 
P : (T,0), we find that 

%(T, 0) = 0, F T% jt k) (T, 0) = 0, % ]M) (T, 0) = 0, • ■ ■ . (A.2) 

Using all these relations regarding the connection coefficients at Po, the re- 
quired derivatives of the metric tensor at Po can be calculated. The results 
can be expressed as 

F g 00 = -1 - F R l0j oX l Xi - l F H iajQh X i XiX k + • • • , (A.3) 

F 9oi = ~ F R Qj ikX j X k + ^ F TZ Qjki£ X^X k X t + ■■■ , (A.4) 

F 9ij = $ij ~ ~^ F RikjeX k X e + - F TZie r jkX k X e X r + • ■ ■ , (A. 5) 

where TZ^p^ is given in terms of covariant derivatives of the Riemann tensor 
by 

T^pLUpcru) "^\Rpupa;ui Rppu)u;cr) ■ (A .6) 

It is simple to show that the tensor 7£ Mi , p(ja , satisfies the following relations 

T^pVpOLQ fipVUXJ J TZpVpOLQ T^fjLpwVO Rp<jp(u;Ul) ■ V^"^) 

The first covariant derivatives of the Riemann tensor are not completely 
specified by the tensor 1Zp VpauJ . Moreover, the next terms in equations (|A.3| )- 
( |A.5j ) contain, among other terms, the tensor 

"TpupcFU-K p,v pau)\-K ~^R^pu)i/Rp,S,ircFi (A .8) 



etc. It follows from these results that the reduced geodesic equation (|2.5| ) in 
Fermi coordinates becomes the tidal equation 

^ + K)X* + ^K) k X^X k + ±K t jkt XiX k X t + • • • = 0, (A.9) 
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where the quantities K^, K^ kl etc., characterize, respectively, the tidal ac- 
celeration terms of the first order, second order, etc. We find from equa- 
tions flOD-(|PD that 



Mjlkj^ 



i£jk^ 



2 ... 

— F R v ^ v v ^ 



as in the generalized Jacobi equation and = 2K^ k s, where 



K' 



ijk 



F n i0j0k + 2[ F n 0jm - -( F n ijm + F n om )\x* - 



+ ( Ttjlikr + g 7tkritj)X X T + -(2 IZojlOk ~ HjOkQl) X % X 



+ -( F 7^nojr — 5 F 7Zek0rj)X l X e X r . 

The Jacobi equation in spaces of arbitrary dimensions was first discussed 
by Levi-Civita 0] and Synge fll], 0, 0- These references contain a detailed 
and explicit exposition of the theorem of Fermi [[13] involving the possibility 
of choosing a system of coordinates in a cylindrical region along an arbitrary 
open curve in spacetime such that all the connection coefficients vanish and 

5V = Viw on tne curve. 

The generalized Jacobi equation was first discussed in jTJ], M and 
developed further in [17|. Later, this generalization was independently redis- 
covered by Ciufolini |TS|, |T9|. Further discussions of the Jacobi equation can 
be found, e.g., in GO, EJ, E3, ffl. 



B Hamiltonian equations of motion 

Consider a coordinate system X a = (T, X 1 ) in the region of interest in space- 
time such that the metric tensor takes the form g^ u = r]^ u + h^ u , where h^ u 
is a small perturbation on the Minkowski background. In particular, we are 
interested in Fermi normal coordinates (J2.ll ). The isoenergetically reduced 



Hamiltonian system is given in these coordinates by 

dX i dH dpi dH 



dT d Vi ' dT dX v ^ R1 ) 
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where 



^T(i±§^) 1/2 (B.2) 
To first order in h^, the Hamiltonian can be expressed as 

H = Tv / lT7(l - \h m - + h 0i Pi , (B.3) 

2 2 1 + p A 

where p 2 = S^PiPj and if is a quadratic form in spatial Fermi coordinates 
X. Using ( |B.3| ), the system of equations ( |B.1| ) can be written as 

dX i _ Pj 1 1 h kt p k pe 



dT 



^^=^~t«> + 2TW )S "- K ' ] + tl "' (B ' 4) 



dT ' 2 v ~ ' r v ~ uu>i T+p 
Letting ^fpi — > pi and defining 



T^V 1 +P 2 (V»+ 1+ 2 )~ h ,iPr ( B - 5 ) 



v 1 + P 



and P* := 5 ij Pj, we find that (1 + p 2 ){l - P 2 ) = 1 and p< = P/Vl - P 2 . In 
terms of P, the Hamiltonian equations of motion ( B.4|) - ([B.5| ) take the form 

d X* 1 1 

= P*(l - 2^00 + ^PfcP) + h oi (B.7) 

1 

2 

^— — d / ^ \ = 1 h _ h 0j p l h M p p (B g) 

dT V^l - P 2 ' 2 ' 2 i* fc * v v 



^Pj + ^h^PkPtP* 



We are interested in the behavior of X % = X l (T) to first order in 5 = 
|X|/g. To this end, we note from equation (B.7|) that P % = V 1 + 0(h^ u ), 
where V™ = dX % jdT. Thus equation ( |B.7| ) can be written to first order in 



P := (1 + \h 00 - + /i Qi + hyVK (B.9) 



25 



Similarly, equation ( B.8 ) can be written as 
+ 



dP % I d~P 1 1 

It + i~^ (p ' dr )pl = 2 Vi + ^ + 2^' iyfey ■ (R10) 



This relation, after multiplying by P{ and summing over i, takes the form 

1 (IP 1, , 1 

dT' 1 ~ 2 



(B.ll) 



1 /'-'•" dT' 2 " ' " UJ "' ' ' 2' 
Substituting equations ( |B.11|) and ( |B.9| ) in equation ( |B.10|) and noting that 

dV i /dT = 0(6), V, = O(5 2 ), 

we finally obtain the desired equation of motion 



dV i 1 

+ {hoojV* + h^y-'V n- - -/loo, 

• (ho;, ~ h tii . ; )\-' + i//,,, - \h lKi )VW h 



0. 



(B.12) 



Taking due account of the approximations under consideration here, we 
note that 



-r° 



dX a dX? 



^HTlT =h ^ V ' + h "^" V 



and 



dX a dX@ 

it,— ^ = r nn + 21^' + r Hk v j 'v k 



jk 



af3 dT dT 



in agreement with equation ( |2.5|) . Moreover, to first order in 5 

hoo,i — —ZRoiQjX 3 , 

2 k 

hoij = —-(Rojik + Rokij)X , 
h^i = -\{Rumk + R mt k)X\ 



in the Fermi system. Substituting these relations in equation ( |B . 1 2| ) and 
using the symmetries of the Riemann tensor, such as Roujk] = 0, we obtain 
equation ( 2.12|) . This is the generalized Jacobi equation in Fermi coordinates; 
for arbitrary coordinates see Appendix C. 
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C Generalized Jacobi equation in arbitrary 
coordinates 

Consider a reference timelike geodesic C as in Appendix A representing the 
worldline of an observer with proper time r. Let = at]' 1 ; then, to first 
order in a/g, the generalized Jacobi equation [14], [L5| takes the form 



n2tM np ncp 



+rK Cpw w — — £ -—J + 2# «/a-FTTS n + ^aTTTS "TtT" ~ U - 



Here 



is the covariant derivative of £ M along the reference worldline. If we now use 
Fermi coordinates and set 

* ~~ (»)' £> r - rf T «' Dr2 _ dT2 A W"' 

we recover equation ( |2.12| ). Here we have used the relation DX'f./Dr = 0; 
i.e. the tetrad is parallel transported along C. 



D Worldlines of observers at rest 

Consider a spacetime metric of the form 

ds 2 = —dt 2 + gij(t, x k )dx l dxK 

If a particle is a rest at a point in the corresponding space, then its worldline 
follows a geodesic. Indeed, the four- velocity of the particle = dx^/dr is 
given by u° = 1 and u % = 0. We must show that 

— + T M u a u p - 
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For fi = 0, it suffices to prove that Tq = 0. Because of the block form of the 
metric, we have g 00 = —1 and g° l = 0; therefore, 

r oo = ~ 2 goo >° = °' 

For n — i, it suffices to show that Tq = 0. In this case we have g j = and 
goo = -1- Hence, 

r oo = ^'(%o,o-goo,i) = 0, 

as required. 
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